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Abstract - We propose a model for the dynamics of a probe embedded in a living cell, where 
both thermal fluctuations and nonequilibrium activity coexist. The model is based on a confining 
harmonic potential describing the elastic cytoskeletal matrix, which undergoes random active hops 
as a result of the nonequilibrium rearrangements within the cell. We describe the probe’s statis¬ 
tics and we bring forth quantities affected by the nonequilibrium activity. We find an excellent 
agreement between the predictions of our model and experimental results for tracers inside living 
cells. Finally, we exploit our model to arrive at quantitative predictions for the parameters char¬ 
acterizing nonequilibrium activity, such as the typical time scale of the activity and the amplitude 
of the active fluctuations. 


Actin filaments are involved in a number of functions 
ncluding cell motility, adhesion, gene expression, and sig- 
lalling. When fueled by ATP supply, Myosin motors 
idvance along these filaments by performing a directed 
itochastic motion. By tracking the trajectory of a micron- 
lize probe embedded within the cytoskeletal network, and 
)y subjecting it to microrheology experiments, one can 
lope to access and understand some of the properties of 
he nonequilibrium activity inside the cytoskeletal net- 
vork. Experiments were first carried out in actin gels 
vithout molecular motors, known as passive gels 
some progress in the experimental field has provided new 
esults for tracers attached to the cortex of living cells , 
md also for in vitro actin gels [^|^. In such gels, called 
ictive gels, the tracer dynamics exhibits large excursions 
corresponding to directed motion events, in addition to 
the thermal fluctuations already observed in passive gels. 
Due to the active processes, the actin network fluctuations 
comprise a strongly nonequilibrium component. Experi¬ 
mentally, the out-of-equilibrium nature of such activity 
has been evidenced by the violation of the fluctuation dis¬ 
sipation theorem (FDT) [9 -11 . To account for nonequilib¬ 
rium activity, a generalization of the FDT has been devel- 
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oped introducing a frequency dependent effective temper- 
This generalization is based on a description 


12-14 


ature 

of tracers dynamics at a mesoscopic scale, which can be de¬ 
scribed using a generalized Langevin equation 15 -17 . At 


a macroscopic scale, the dynamics of acto-myosin networks 
have been described via hydrodynamic treatments 


polymer theory 19 20 
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or 


In what follows, we present results of microrheology ex¬ 
periments in the cytoplasm of living cells, which are char¬ 
acterized by a highly nonequilibrium activity. Along with 
experiments, we propose a model which mixes simple but 
nontrivial rheology with random fluctuations due to active 
processes inside the cell. We carry out a comparison with 
experimental data, which allows us to directly determine 
some microscopic mechanisms that drive active fluctua¬ 
tions inside the cell. We demonstrate that our quantitive 
estimation of the nonequilibrium active features is consis¬ 
tent with different kinds of experimental measurements, 
thus supporting the overall consistency of our model. 


We inject sub-micron colloidal tracers in the cytoplasm 
of living A7 cells, and track a two-dimensional projec¬ 
tion of their fluctuating (3-D) motion with confocal mi- 
We observe some directed motion events in 
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croscopy 

the tracers’ trajectories in addition to the thermal fluctu- 
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Fig. 1: (Top) Typical trajectories of 200 nm PEG coated beads 
in A7 cells under three conditions: control, 10 pM blebbistatin 
treatment, ATP depletion. Trajectory length is about 2 min. 
(Bottom) Elastic storage modulus G' (+) and loss modulus 
G” (o) from active microrheology experiments in untreated 
(red), blebbistatin treated (orange), and ATP depleted (blue) 
A7 cells. 


ations of small amplitudes (Fig. [^, as already reported 
in synthetic active gels . To investigate the role of bio¬ 
logical activity in the intracellular mechanics, we subject 
cells to two treatments. We inhibit Myosin II motors by 
adding 10 /iM of blebbistatin to the culture medium, and 
we deplete cells of ATP through addition of 2 mM sodium 
azide and 10 mM of 2-deoxyglucose. We extract the one¬ 
dimensional mean square displacement (MSD) from the 
spontaneous motion of tracers for different radius sizes 
a = {50,100, 250} nm. We present the MSD multiplied 
by a for the control, blebbistatin and ATP depleted con¬ 
ditions in Fig.j^a), showing that the MSD scales like 1/a. 
The small time MSD is constant in the three conditions, 
while the large time behavior is diffusive, apart for ATP 
depleted cells, where it remains almost constant. Since the 
time evolution of the MSD is qualitatively similar for trac¬ 
ers of different sizes, we deduce that the tracers are bigger 
than the mesh size of the cytoskeletal network, thus allow¬ 
ing us to consider that they evolve in a continuous medium 
in first approximation. 

We measure the mechanical properties of the cytoplasm 
via active microrheology method by using optical tweez¬ 
ers 22 . We impose a sinusoidal oscillation on a particle 
with diameter 0.5 fim within the cytoplasm. From the 
resultant displacement of the bead, we extract the com¬ 
plex modulus G* = l/(67rax), where x is the Fourier re¬ 
sponse function. It reveals that it weakly depends on fre¬ 
quency, and that the elastic contribution is significantly 
larger than the dissipative one (Fig.[^, in agreement with 


previous results 22 23 . Moreover, we do not observe 


a significant change in the cytoplasmic mechanical prop¬ 
erty according to active processes. The cytoplasm is still 
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Fig. 2: (a) Time evolution of the one-dimensional mean square 
displacement scaled with the tracer radius a = 50 (-h), 100 
(o) and 250 nm (•) for control (red), blebbistatin treated (or¬ 
ange), and ATP depleted (blue) cells, and the corresponding 
best fitting curves (Eq. @): solid, dashed, and dot dashed line, 
respectively. (Inset) FDT-ratio as a function of frequency. It 
eqnals 1 in ATP depleted cells as for an equilibrinm system, 
and it deviates from it in the two other conditions at small 
frequency showing that nonequilibrium processes drive the dy¬ 
namics in this regime, (b) Time evolntion of the MSD scaled 
with G' measnred with tracers of radins a = 100 nm in control 
cells. The G' value increases with the percentage of PEG in¬ 
troduced in the cell: 0% (yellow •), 3% (light green o), and 6% 
(dark green -I-). The best fit curves are in solid, dashed, and 
dot dashed lines, respectively. The short time scale plateau 
scales like 1/G'. (c) Time evolution of the MSD times G'^. 
The large time diffusive part scales as IjG'^. 


mainly elastic in blebbistatin treated and ATP depleted 
cells, with a storage modulus being twice as small as in 
untreated cells where it equals approximately 2 Pa. 

To quantify departure from equilibrium, we extract the 
FDT-ratio which compares the active microrheology mea- 


p-2 






















Activity driven fluctuations in living cells 


surement with the random intracellular motion visual¬ 
ized by tracer particles [9 11 . It is defined in terms of 
the position power spectrum C and the imaginary part 
of the Fourier response function y" as FDT-ratio(a;) = 
—ojC{oj)/ [2y"(a;)fcB'r], where T is the bath temperature. 
It equals 1 for an equilibrium system, and deviates from it 
otherwise. The control and blebbistatin treated cells are 
out-of-equilibrium, whereas the effect of the nonequilib¬ 
rium processes are negligible in ATP depleted cells (Inset 
in Fig.|^a)). This supports that the nonequilibrium pro¬ 
cesses hibernate in the latter as long as no ATP supply 
is provided, suggesting that there is an equilibrium refer¬ 
ence state where the tracer particle is trapped in an elastic 
cytoskeletal network. Given we can not rely on equilib¬ 
rium physics to describe the tracer’s dynamics in the two 
other conditions, we offer a new model to characterize its 
nonequilibrium properties. 

We vary experimentally the elastic modulus G' by 
adding various amount of 300-Dalton polyethylene glycol 
(PEG) into the cell culture medium Q This results in an 
osmotic compression on the cell, so that G' increases with 
the amount of PEG applied [^. We report in Figs.[^(b)- 
(c) the MSD data multiplied by G' and G'^ for different 
values of G'. It appears the value of the small time plateau 
scales as 1/G' while the large time diffusion constant scales 
as l/G' 2 . 

The cytoskeleton acts as a thermostat for the tracer par¬ 
ticle. Provided that inertial effects are negligible in the 
intracellular environment, we model the dynamics of the 
tracer’s position r by means of an overdamped Langevin 
equation. We use a harmonic approximation to account 
for the interaction of the tracer with the surrounding net¬ 
work. The main new ingredient of our model lies in ex¬ 
pressing the effect of nonequilibrium activity. We postu¬ 
late that the underlying action of the active processes in¬ 
duces local rearrangements of the network, resulting in an 
active force applied on the tracers. As an example of such 
nonequilibrium processes, the activity of Myosin II motors 
can slide cytoskeletal filaments past each others leading to 
a local deformation of the network . To account for the 
directed motion events observed in our experimental tra¬ 
jectories, we consider that the active force proceeds by 
a sequence of rapid ballistic jumps followed by quiescent 
periods. It remains constant during intervals of average 
quiescence time tq, when the tracer is only subjected to 
thermal fluctuations, and it varies during a persistence 
time of order r by a quantity Ta = /n, where n is a ran¬ 
dom direction in the three-dimensional space. We assume 
that the persistence and quiescence times are exponen¬ 
tially distributed variables as observed in synthetic active 
gels [7|[^[^ , and that they do not depend on the network 
and tracer properties. Putting these ingredients together, 
we arrive at the equation for cc, the one-dimensional pro- 


^ After the stress, cells are allowed to equilibrate for 10 min at 
37° C and 5% CO 2 , before we perform the imaging or optical-tweezer 
measurement. The cell size and mechanics equilibrate in 2 min after 
adding PEG based on our imaging and previous studies [24|. 



Fig. 3: Typical realization of (a) the active force /a, and 
(b) the corresponding active bursts Va- /a is constant over 
a quiescence time of typical value tq, and varies linearly with 
a slope uniformly distributed in [—/, /] during a persistence 
time of order r. Va is proportional to the time derivative of 
/a. (c) Schematic representation of the energetic landscape re¬ 
arrangement due to nonequilibrium activity and its modeling 
using the active burst applied on the local minimum. We depict 
the network potential in black solid line, the harmonic approx¬ 
imation in dashed red line, and the tracer particle in Hlled blue 
circle. Nonequilibrium activity leads to a displacement vr of 
the potential, resulting in an energy gain E ~ k{vT)^ for the 
tracer. 


jection of r 


dx , A r 

1-^ = -kx + t, + Ja, 


( 1 ) 


where ^ is a zero mean Gaussian white noise with correla¬ 
tions = 2jkBTS{t—t'), and /a is a random force 

with typical realization described in Fig. [^a). The spring 
constant of the surrounding network is k, and 7 is the 
friction coefficient of the environment. Our model is asso¬ 
ciated with a Fourier response function y = l/(k + iuij), 
from which we deduce that complex modulus is of the form 
G* = l/( 67 rax) = k/{6Tra) +iu!ri, where 77 is the viscosity 


of the fluid surrounding the tracer 17 2^. We neglect the 


weak frequency dependence of the real part G' as deter¬ 
mined from active microrheology measurements, so that 
the spring constant is directly given by fc = dnaG', as al¬ 
ready reported in other complex fluids with similar elastic 
behavior 26 . Stokes’ law ensures that 7 is independent 


of G', and 7 oc a. 

To illustrate our model with an immediate physical pic¬ 
ture, we introduce the variable Tq = tA/k which we regard 
as the position of the local minimum of the potential in 
which the tracer is trapped. The local rearrangements 
of the network due to nonequilibrium activity result in a 
shift of the local minimum the tracer sits in. Thus, this 
position has a dynamics of its own given by a random 
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active burst Va in which a burst r;n is felt during the per¬ 
sistence time, while it equals zero during the quiescence 
time (Fig. [^b)). The active force projection is simply re¬ 
lated to the active burst projection as dfj^/dt = kvj^. We 
assume that the typical variation / of the active force is 
independent of the network properties, whereas the active 
burst amplitude v = //(kr) depends on the properties of 
the cytoskeletal network via k. 

From the Fourier transform of Eq. 0. we compute the 
position autocorrelation function C{t) = {x{t)x{0)), and 
then deduce the one-dimensional MSD as (Aa;^) (t) = 
2((7(0) — C{t)). We denote the thermal contribution to 
the MSD by (Ax^), and the MSD when the particle 
is only subjected to motor activity by (Ax^), so that: 
(Ax^) = (Ax^) -I- (Axa). The thermal MSD is the same 
as for the Ornstein-Uhlenbeck process, and we compute 
the active contribution in terms of the parameters charac¬ 
terizing the active force: 


(Ax^) (t) 
(Axl) (t) 


2k^T 
k 

2kjiTj^/k 

1 - {t/t.Y 


(i 


( 2 a) 




-be' 


-t/r^ 


- 1 


( 2 b) 


where = y/fc is a microscopic relaxation time scale. In 
the passive case, i.e. when Ta = 0 , it saturates to the value 
2k^T/k within a time Tr as predicted by the equipartition 
theorem, meaning that the tracer is confined in the cy- 
toskeleton. The active force represents the random fluctu¬ 
ations of the cytoskeletal network induced by the nonequi¬ 
librium activity. With such a force, the MSD exhibits a 
plateau at the equilibrium value corresponding to a tran¬ 
sient elastic confinement at times Tj- t t, and then 
has a diffusion-like growth on longer times with coefficient 
2k-B,Tj^/^. Provided that k oc G', it follows that the equilib¬ 
rium plateau scales like IjG', as we observe experimentally 
(Fig. [^b)). The energy scale k^Tj^ = 7 (xr)^/[ 3 (T -b tq)] 
defines an active temperature, which is related to the am¬ 
plitude of the active fluctuations as defined by the ac¬ 
tive burst correlations {v/ir^). 
The independence of / and r with respect to G' yields 
V oc 1/G', from which we deduce that the large time dif¬ 
fusion coefficient scales as XjG''^^ in agreement with our 
measurements (Fig. 0 c)). 

On the basis of our phenomenological picture where the 
nonequilibrium dynamics is driven by an active remod¬ 
elling of the cytoskeletal network, we propose a physical 
argument for the scaling of the MSD with the tracers’ size 
a presented in Fig. 0a). As presented above, we first as¬ 
sume that k and 7 scales like a. Within our model, the 
active burst represents the activity-driven network defor¬ 
mation and reorganization, which result in a change of the 
tracer’s local energetic landscape. During a burst event, 
the local minimum is shifted by a random amount. Re¬ 
garding this event as instantaneous, the tracer finds itself 


at a distance of order vt from the new local minimum 
position after each burst. It follows that the typical en¬ 
ergy provided by nonequilibrium activity to the particle is 
E ~ k{vT)'^, as depicted in Fig. 0c). We assume that it 
does not depend on the particle properties, just as t and 
To, thus being independent of the tracer’s typical size a. 
Since fc oc a, we deduce v cx implying that Ta is 

independent of a. Finally, the relaxation time Tr is also 
independent of a, leading to a scaling of the MSD like 1/a 
which agrees with our observation. 

We use our analytic expression to fit the MSD data mul¬ 
tiplied by a for the three conditions described above. We 
assume the viscosity of the fluid surrounding the tracer is 
the cytoplasm viscosity rj ^ 10“^ Pa • s 


27 


and we de¬ 
duce the damping coefficient from Stokes’ law: 7 = Girarj. 
We estimate the k value from the small time plateau. 
The only remaining parameters are the ones character¬ 
izing nonequilibrium activity: Tj^/T = {2.8,0.9} x 10“^, 
r = {0.16±0.03, 0.39±0.09} s, in control and blebbistatin 
treated cells, respectively. The estimation error made on 
Ta/T is of the order of 1 % in control, and 0 . 1 % in bleb¬ 
bistatin treated cells. 

The amplitude of the active fluctuations is smaller in 
blebbistatin treated cells, meaning that the inhibition of 
Myosin II motors reduces the proportion of nonequilib¬ 
rium fluctuations with respect to the thermal ones as ex¬ 
pected. Other nonequilibrium processes drive the out-of- 
equilibrium dynamics in this condition. The typical time 
scale r of the persistent motion events is enhanced in bleb¬ 
bistatin treated cells. Assuming that each active burst 
persists until the stress that accumulates in the network 
causes the network to locally fail, weaker motors due to 
the addition of blebistatin will contract for a longer dura¬ 
tion until such a critical stress builds up. Provided that 
l/r is the typical frequency below which the nonequilib¬ 
rium processes affect the dynamics, this supports that the 
active fluctuations take over the thermal ones at larger 
frequencies in the control cells compared with the blebbis¬ 
tatin treated ones. Notice that Ta represents the ability 
of the tracer to diffuse on long times, and T quantifies 
here only the motion of the bead at short times when it is 
trapped within the elastic cytoskeletal network. The fact 
that we find Ta small compared to T does not mean that 
the active processes are negligible, as they control entirely 
the long-time and long-distance diffusion of the tracer. In 
the absence of activity, the tracer does not diffuse at all 
and remains trapped in the elastic network. 

To characterize the properties of the active force, we 
focus on the power spectrum of the stress fluctuations, 
i. e. the Fourier transform of the time correlation function 
(/A(0)/A(t)) 0 We extract the power spectrum 

of the overall force /a + C as the power spectrum of the 
position times ( 67 ra|G*|)^ 15 . Provided that the ATP 
depleted condition is in an equilibrium state, the active 
force /a is negligible in these cells and the overall force 
reduces to the thus providing a direct measurement of 
the thermal force spectrum. We remove this equilibrium 
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contribution to the overall spectrum to deduce the active 
force spectrum in the two other conditions. We observe 
a l/w^ behavior at low frequency as already accounted 
for on general grounds [^[^[^[^, and the large fre¬ 
quency curvature hints a crossover to another power law 
(Figs. I^a)-(b)). Our analytic prediction for the active 
force spectrum reads 




1 2j keTj^ 
(wTr)^ 1 -I- (wr)^ 


(3) 


It combines properties of the network and parameters 
characterizing the active force, since the effect of nonequi¬ 
librium activity on the tracer is mediated by the network 
within our model. We recover the divergence as XjuP' at 
low frequency, and we predict a power law behavior l/w^ 
at high frequency, the crossover between the two regimes 
appearing at 1/r. We compare our prediction with the 
experimental data by using the best fit parameters esti¬ 
mated from the MSD data. Without any free parameter, 
we reproduce the measured spectra (Figs. |^a)-(b)). This 
result is a strong support for our model, in which not 
only quantifies the long time diffusion coefficient of the 
tracers, it is also related to the typical amplitude of the 
fluctuations generated by the nonequilibrium active force. 
The study of the high frequency spectrum calls for new 
experiments as it would confirm the validity of our phe¬ 
nomenological picture. 

To study in more details the properties of the active 
force, we analyze the probability distribution function of 
the tracer displacement (DPDF). It exhibits a Gaussian 
behavior at short and long times. In the intermediate 
regime, we observe a central Gaussian part which matches 
our equilibrium prediction in the absence of activity, and 
exponential tails accounting for directed motion events 
consistent with previous observations in synthetic active 
gels [^. Within our model, the non-Gaussian behavior 
of the DPDF is a direct and unique consequence of the 
non-Gaussianity of the active force. We ran numerical 
simulations of the dynamics in Eq. Q to reproduce the 
time evolution of the DPDF. We set the different parame¬ 
ter values to the one estimated previously, letting us with 
only one free parameter: the average quiescence time tq. 
It quantifies the average time between two successive di¬ 
rected motion events, thus controlling the relative impor¬ 
tance of the exponential tails with respect to the Gaussian 
central part. We adjust this parameter by matching the 
exponential tails observed at different times. 

With a fixed tq value, we manage to reproduce the evo¬ 
lution in time of the whole experimental DPDF. This 
shows that the specific form we choose for the active 
process is sufficient to reproduce not only the MSD and 
force spectrum data, but also to account quantitatively for 
the dynamic non-Gaussian properties of the distribution 
(Figs. i c)-(d)). We estimate tq = {2.5, 2.8} s in control 
and blebbistatin treated cells, respectively. The extracted 
values are very similar for the two conditions, showing 




Fig. 4: (Top) Active force spectrum Sa as a function of fre¬ 
quency measured with tracers of radius a = 250 nm in (a) con¬ 
trol , and (b) blebbistatin treated cells. The experimental data 
are in black • , and the dashed lines correspond to Eq. with 
the parameter values deduced from the best fit of the MSD 
data. (Bottom) Probability distribution function of the tracer 
displacement (DPDF) at two different times: 1 (o), and 5 s (•). 
The DPDF is measured with tracers of radius a — 250 nm in 
(c) control, and (b) blebbistatin treated cells. We present the 
corresponding results from numerical simulations of Eq. 0 in 
solid and dashed lines, respectively. The blue dot dashed line 
is the corresponding equilibrium Gaussian. The parameter val¬ 
ues are the same for the two lag times: (b) {Ta/T, r, tq, k, 7 } = 
{2.8 X 10"®,0.16s,2.5s,8.5pN/pm,4.7 x 10"® pN • s/pm}, 
(c) {Ta/T, t. To, fc, 7 } = {9 X 10“'‘, 0.39 s, 2.8 s, 8.2 pN/pm, 4.7x 
10“® pN • s/pm}. 


that the addition of blebbistatin does not affect the typical 
time over which the tracers are only subjected to thermal 
fluctuations. It suggests that this time scale is related to 
the recovery of the network following a large reorganiza¬ 
tion, thus being barely independent of the activity of the 
nonequilibrium processes. Notice that the corresponding 
duty ratio Pon = t/{t + tq) is smaller in control than 
in the blebbistatin treated cells: Pon = {6,15}%, respec¬ 
tively. It is a quantitative evidence that the exponential 
tails are more pronounced in the control condition, namely 
the proportion of directed motion events is increased. We 
deduce the value of the typical active burst amplitude: 
V = {0.86,0.22} pm/s in control and blebbistatin treated 
cells for a = 250 nm, which are compatible with velocity 
scales observed in [^ . 

Microrheology methods have become a standard tech¬ 
nique to explore cellular activity in living organisms 30 


P-5 









E. Fodor et al. 


In this work, we introduce a new model for characterizing 
the motion of a tracer in a living cell. This model explic¬ 
itly accounts for the elastic behavior of the cytoskeletal 
network and successfully combines it with a description 
of the cellular active force —a well defined non-Gaussian 
colored process. By analyzing the MSD data, we quantify 
two essential features of this force: its strength, and the 
typical time scale over which it is felt. Our model goes be¬ 
yond previous modeling which treated the nonequilibrium 
activity as a random noise with unprescribed character- 
In a previous work, activity was modeled as 


15 


istics 

a trichotomous noise acting directly on the particle 12 


whereas such activity is mediated by the surrounding net¬ 
work within our new proposal. The present model com¬ 
bines the short time confined behavior with a long time 
free diffusion which is driven by the active force, and re¬ 
covers all the main experimental results. The model ap¬ 
plies as long as we are in the regime of simple viscoelastic 
behavior. Dressing our model with a more realistic rhe¬ 
ology, e.g. with a power law behavior for the complex 
modulus, usually observed in cell rheology 31 , is concep¬ 


tually straightforward as a future elaboration of the model. 
Further generalization of our model could be used to de¬ 
scribe active fluctuations in other non-equilibrium (living 
or mechanically driven) systems that exhibits similar be¬ 
havior [7 32 

We acknowledge several useful discussions with Julien 
Tailleur and Frangois Gallet. N.S.G. gratefully acknowl¬ 
edges funding from the Israel Science Foundation (grant 
no. 580/12). 
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